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SUMMARY 


A theoretical analysis is made of a square simply supported plate 
loaded in edge compression and subject to creep. The plate Is assumed 
to be made of a material which obeys a- nonlinear creep lav. An approxi- 
mate solution is carried out with the use of a previously published vari- 
ational theorem for creep. The theory does not yield a finite collapse 
time but does indicate how lateral deflections and unit shortening due 
to creep might be calculated. The results also show that creep can 
cause significant redistribution of the middle -surface stresses in a 
plate . 


INTRODUCTION 


At elevated temperatures, aircraft structural materials may exhibit 
the phenomenon of creep (deformation with time at constant load). The 
need for aircraft to operate in an elevated temperature environment, 
therefore, has Increased the importance of considering creep in their 
structural design. As a result, interest has been stimulated in the 
creep behavior of such structural components as plates subjected to 
in-plane edge loads. 

Some experimental and theoretical work, has been published on rec- 
tangular plates subject to creep. In reference 1 results of creep tests 
and an empirical method for predicting collapse times are presented for 
plates loaded in compression on two opposite edges and with the remaining 
edges unloaded and supported in V-groove fixtures. Other approximate 
methods for handling plates having various types of edge support are 
suggested in reference d; however, experimental verification for these 
methods Is quite limited. 

In reference 3 an analysis based on small-deflection theory is made 
of the creep deflection of a simply supported plate composed of a linear 
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viscoelastic material - that is, a material in which the stress and 
strain and their appropriate time derivatives are related In a linear 
fashion. Unfortunately, most structural materials do not exhibit the 
properties of linear viscoelasticity. The creep relation Is usually 
nonlinear; for example, it is sometimes expressed by a power relation- 
ship between strain rate and stress. 

A variational theorem given in reference 4 provides a way of over- 
coming some of the difficulties caused by the nonlinear creep law. This 
theorem facilitates approximate analyses for plates of no nline ar mate- 
rials by the direct methods of the calculus of variations. One such 
analysis utilizing this variational theorem is presented in reference 5 
for a sandwich plate. A solution is given in the present paper for a 
square simply supported solid plate loaded in compression along two 
opposite edges. In this investigation all edges of the plate are 
required to remain straight and in the original plane as the plate 
deflects. 


SYMBOIS 



F 

f 

h 

I 


width or length of square plate 
effective width of -plate 
Young 1 s modulus 
unit shortening of plate 


unit shortening at elastic buckling load, 



function of Jg and t (see eq. (10)) 
unit displacement of plate side edges 
thickness of plate 

denotes various integrals defined in equa- 
tions ( 25 ) 
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J 2 


k 

n 

P 



■^cr " CT crbk 
S 

&u 

S a 

8 ij 

*i 

*1 

t 


second invariant of stress deviation tensor, 

| s iJ s ij 

constant in creep lav (see eq. (ll)) 

exponent in creep lav (see eq. (ll)) 

total compressive load on plate, positive in 
compression 


area of a surface 

surface where displacement rates are prescribed 
surface where stress rates are prescribed 

stress deviation tensor, cr-y - ^ 

surface traction 
prescribed surface traction 
time 



u,v,v displacements of the middle surface in x-, y-, 

and z-directions, respectively 

displacement vector 

u 10* u ll' v 01> v ll^ W 11 arbitrary coefficients in equations ( 19 ) 


V 


volume 


W = 


V ll b 

h 


Wi 


v iii b 

h 


amplitude coefficient of initial imperfection 
of plate (see eq. (A4)) 


w lli 
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x,y,z 

x i 

^xy-’^xz-’^yz 

8 ij 

e U 

€ ij 

II 

e i«J 

e x> € y* e z 

?xy 
€ xy “ — 





A 

H 

n 

% 


plate coordinates 

plate coordinates in tensor notation 

shear strains in xy-, xz-, and yz -direct ions, 
respectively 

Krone cker delta 
strain tensor 

noncreep portion of strain tensor 
creep strain tensor 

strains in x-, y-, and z-directions, respectively 




•rfy 3CV* Zt 

dimensionless coordinates; 
respectively 

constant in creep lav (see eq. (ll)) 

Poisson's ratio 

quantity to he varied in a variational theorem 

quantity to he varied in Reissner's variational 
theorem 

elastic buckling stress for square plate, 

R /jth\2 1 

3(1 - ^ 2 ) 
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J i<5 


stress tensor 


o x , cry- 


stress in x- and y-direetions, respectively, 
positive in tension 


°xOOj ^xlO; °x01^ ff 3dlJ °xb>GyOO> 

°ylO-’ °y01> °yll' a yb jT xyll^ T xyb 


l arbitrary coefficients in equations (20) 


time parameter. 


^ Ea cr 2n t l£ 


®,f,e 

Bar denotes 

Dot denotes differentiation with respect to time. 

Repeated_letter_subscrl.pt denotes s umma tion over 1, 2, 3j for 
instance, u-jT^ = upT^ + ugTg + u^Tj. 

Ccamna preceding a subscript denotes partial differentiation with 

• 8uv 

respect to x with that subscript; for instance, uv j_ = — — . 

OX-j^ 


shear stress in planes parallel to xy -plane 

quantities defined in equations (22) 
a dimensionless quantity. 


ANAETSIS 


Statement of Problem 

The plate considered in this investigation is illustrated in fig- 
ure 1. It is a square flat plate compressed in its plane along the two 
edges parallel to the y-axis. All edges are assumed to be simply sup- 
ported and to remain in the original plane of the plate. The loaded 
edges are assumed to be constrained to remain straight and to be free 
of tangential stresses; that is, the plate is loaded by means of a pair 
of rigid frictionless loading platens. The unloaded edges also are 
assumed to be constrained to remain straight and free of tangential 
stresses and are free to translate in the plane of the plate. The load 
application is assumed to be sufficiently rapid so that the creep which 
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occurs during loading is negligible but sufficiently slaw so that dynamic 
effects also are negligible. After the total load is applied it is 
assumed to remain constant with time. The history of the lateral deflec- 
tion, unit shortening, and middle-surface stresses are desired subsequent 
to load application. 

The boundary conditions, which are considered to hold for all time, 
can be written in mathematical form as follows: 


For constant displacement, 



For zero deflection, 


v^±|,yj = v^x,±|j -0 

For zero moment. 



For loaded edges. 



( 1 ) 


( 2 ) 


( 3 ) 


( 4 ) 


For unloaded edges. 


r h /2 pV 2 
o -h/2 J -b/2 ' 



dx dz 


= 0 


( 3 ) 


For shear stress, 
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Variational Theorem 

In this analysis the variational theorem presented in reference 4 
is utilized. Thi s theorem differs from the usual variational theorems 
for potential energy and complementary energy in that instead of varying 
stress and strain the stress rates and strain rates are varied. In this 
respect the theorem is similar to those presented in reference 6. Fur- 
thermore, in the variational theorem of reference 4, stress rates nnH 
displacement rates can be approximated independently) in this respect 
the theorem is similar to that of Eeissner in reference 7* The theorem 
of reference 4 is stated here in general form and subsequently reduced 
to a form appropriate for the present problem by incorporating strain- 
displacement relations and a creep law. 

General statement of theorem .- The variational theorem of refer- 
ence 4 is stated as follows : If the stresses and strains are known 

throughout a body at a given instant of time, then the stress rates 
and strain rates existing at that instant are given by 811 = 0 where 


n ■ X + 2 ^«k,J ff i J " + 


dV - 



(h - 38 


(7) 


In the first term of the volume integral, is expressed in terms 

of displacements and displacement rates. In the third term of the 
volume integral, is expressed In terms of stresses and stress 

rates, and €j_j is expressed in terms of stresses and time and is 

independent of stress rate. The variations are to be carried out only 
with respect to time derivatives of quantities . 

Strain-displacement relations . - It Is assumed that lines which are 
originally normal to the middle surface remain normal during the deflec- 
tion of the plate. As a consequence, the strain-displacement relations 
can /be written as follows : 
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V 


8v + lj^dw'' 

8y 



- z 


S 2 w 

dy 2 


2e xy y 


xy 


= £n + 8v t 8w 8w 
8y 5x Sx Sy 


2z 


*2 
O V 

Six Sy 


e z ~ 7xz ?yz 


= 0 


(8b) 

(8c) 

(8a) 


Equations (8) are the strain-displacement relations consistent with the 
Von Kn.rman equations for large deflections of plates. The time deriva- 
tives of equations (8) are 


?xy 



6 X = 

! + 

Sw 

Sx 

Sw S 2 w 


(9a) 


*y = 


Sw 

8y 

Sw S 2 w 

Sy 2 


(913 ) 

_ 8u 


, Sw 

Sw 

Sw Sw p 

-s2. 

0 V 

(9c) 

8y 

+ — 
Sx 

Sx 

dy 

T n l—Li 

Sx Sy 

Six Sy 



k z 

?xz 

11 

0 


(9d) 


Creep lav .- In the formulation of the variational theorem, It is 
assumed that the total strain can be separated into an elastic-plastic 
part and a creep part . Plastic strains are neglected in the 

present investigation; thus, Is given in terms of stress by Hooke' 

law. The functional relation between the creep rate tensor and the 
stress tensor and time may be assumed as fallows: 
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where s-y is the stress deviation tensor and is related to the stress 
tensor by 


s ij " a ' 


ij “ 3 °kk 5 i«3 


Equation (10) reflects the fact that creep strains appear to be inde- 
pendent of hydrostatic stress. The form of F is taken to be 


F = AJ 2 n tt k " 1 


( 11 ) 


where A, n, and k are material constants . The symbol J 2 repre- 
sents the second invariant of the stress deviation tensor and is given 
by J 2 = s^js^j. The power function of <J 2 is suggested in refer- 
ences 8 and 9 where some experimental support is presented for it. The 
power function of time provides an approximate way to account for primary 
creep. 


Equation (10), of course, applies to three-dimensional creep prob- 
lems. A two-dimensional form applicable to plates can be obtained when 


= a z = 0 


and 


°i 3 = T xz “ a 23 = T yz “ 0 


and it is written as follows in engineering notation: 


e ll “ e x 


_ ( g x 2 + °y 2 ~ + 5 T xy 2 ) 


■(20x “ Cfy)kt 


k -1 


(g x 2 + Oy 2 - J X Jy + ^Txy 2 ) ^ (12) 


. if .it A 

*22 = ®y = 7 

•" A (°x 2 + °y 2 “ °x°y + ^ T xy 2 ) 

2£ 12 = V - 5 T 

5 


■(^Txy)^ 


k -1 
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In the one-dimensional case the additional stresses = 0 y and 

r j _2 = T xy W set equal to zero, and the following equation results: 

4 = f * ktk_1 ( 13 ) 

5 3 n 

Form of theorem appropriate for present problem .- With the use of 
the strain-displacement relations and the creep law presented in the 
previous sections, the variational theorem can he reduced to a form 
appropriate for the problem under consideration. When equations ( 9 ) are 
used, the first term in the volume integral of the equation for II 
(eq. (7)) becomes 


eijtfij = r- 


9u 9w 9w 

8x 9 k 9x 


9v 9w 9w 9 2 w\. 

+ Z Idy + 

dy dy $y N 2/ 


chi , 9v , 9w 9w , 9w 9w p _9jw_\ • 

9y 8x 8x 9y 9x 9y 9x 9yj ^ 


Notice that to obtain equation (l^+) the fact that T^y = is used. 

The second term of the volume integral in equation ( 7 ) yields the fol- 
lowing expression when certain terms are neglected to be consistent 
with equations (8) : 


1 |/9w-\ 2 


2 J a ij = 2 Is r x + r 7 + sz ^ Tx y 


9w 9w 


The third term of the volume integral becomes: 


1 • ' u _ 1 /1 2 , i 2^ n i : , I + 11 • 

2 4 y 2E\ X y / E xy E 3 


The final term in the volume integral, the creep term, becomes: 

^I^ij = ^ *^2^) [K - °y )5 ^ + (2o 3r “ a x) b y + ^xy*xy] (17) 



MCA TN 4398 


11 


when equation (10) is utilized to express in terms of stress 

and time. 

The two surface integrals in the equation for IT (eq. ( 7 )) remain 
to be considered. Displacement rates are not prescribed anywhere on 
the boundaries of the plate. The surface integral over Su, therefore, 
may be discarded. The other surface integral, the one over Sq, must 
be modified to account for the loading condition and the requirement of 
straight edges. The correct modification is to replace the surface 
integral over S 0 by the term Peb where P represents the total load 
on the plate and e represents the plate unit shortening. 

When equations (l4) to ( 17 ) are substituted into the equation for II 
(eq. ( 7 )) and the correct modification of the surface integral over S 0 

is used, the following expression appropriate for the problem under 
consideration results: 


n = 




dv [ dw 5w ; dvr dw 
Sx cbc by 8x 8y 


dw dw 
8x &x 



- 2z 



+ 


8v t 8w 8w 
Sy 8y 9y 






f* t™. - A-A 2 + a 2 ) r 

8x by ** 2EV y / E xy E 

1 


Sy 

i f|^2o- x - ffy) She + (20y - cr x )c y + Gt^t 


dx dy dz - Peb 


( 18 ) 


Method of Solution 


A Rayleigh-Ritz type procedure can be used in conjunction with 
equation (l8) to obtain an approximate solution to this problem. Reason- 
able assumptions are made for the stresses, displacements, stress rates, 
and displacement rates in the plate. These assumed distributions contain 
certain arbitrary parameters and the time rates of change of these arbi- 
trary parameters. When the assumed distributions are substituted into 
the expression for H, integrations over the volume of the plate can be 
carried out. Then a system of equations is obtained by equating to zero 
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the variation of the resulting expression for II with respect to varia- 
tions in the time rates of change of the arbitrary parameters. These 
equations are, in general, nonlinear ordinary differential equations of 
the first order with time as the independent variable. Numerical methods 
are available for the solution of such systems when appropriate initial 
conditions are known. 


Derivation of Differential Equations 

As siimed stress rates and displacement rates .- The boundary condi- 
tions considered in this report are the same as those of reference 10. 
Therefore, the displacement functions assumed in reference 10 are used 
as a guide in choosing the displacement-rate and stress-rate distribu- 
tions. The middle-surface displacement rates are assumed to be the 
time derivatives of the displacements of reference 10: 


u = -ex + b Uj_o + U -; -| cos 




* / 


v = fy + b 


(*01 + 4 


T 11 ooa 


2itx 


* / 


I Bin 


(19) 


• , • Jtx Jty 

w = bw-n cos — cos — 

J - L b b 


In choosing appropriate stress-rate distributions, some assumption 
must be made regarding the variation of stress through the thickness of 
the plate. As pointed out in reference k, this variation Is not neces- 
sarily linear even though the strain variation is assumed linear through 
the thickness. Calculations of creep of columns, however, suggest t ha t 
in certain cases the nonlinearity of the stress distribution may not 
seriously Influence the lateral deflection and unit shortening. On this 
basis. It is assumed that a linear variation of stress rate with z Is 
a sufficiently accurate first approximation for plates, and the stress 
rates are taken as follows: 


°x “ ®x00 + ®xl0 006 ~r— + ®x01 cos + 0*11 ° os -TT- 008 yy + oib 00 s ry cos 


• • • 2jcx • 2jcy 

°y = °y00 + °y!0 cos — + cya. cos — £ + 


at* 

b 


gaty 

~b~ 


• 

2z 

JtX 

Jty 

<*rb 

b 

COB — 
D 

cos 

D 


2z 

3fX 

xy 

Cyb 

b 

COS 

■h 

cos — 
■»* 


1 2*x 2*y • 2z . joc . «y 

T— y ^3Cvll Sin ' + Sin Sin ■" 

^ * yj - L b b xyD h. b b 


( 20 ) 
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The displacement rates (eqs. (19)) and "the stress rates (eqs. (20)) 
satisfy the time derivatives of the boundary conditions given by equa- 
tions (l) to ( 6 ). It is in these assumptions for the stress rates and 
displacement rates that this investigation differs from that of refer- 
ence 5 . In reference 5 the simplified assumptions are made that the' 
stresses are related to the lateral deflection in the same manner as 
that of the elastic buckling of a plate under compression in one direc- 
tion. This procedure allows all, the unknown coefficients to be written 
in terms of w-^, the lateral deflection coefficient. The variation 

of n, then, leads to a single differential equation in w-^ instead 
of the system of simultaneous differential equations in Wj ,u-[ n , u-, , . . 
a xOO> a xlO> ‘ * * ■wkick result from the present analysis. 

Variation of II. - When equations (19) and (20) are substituted into 
equation (l 8 ) and certain of the integrations are carried out, the fol- 
lowing expression results: 


. A 


.. • “ 11 . Tu. • “ 11 + 

-=°iOO + f ®y00 + *l«UO®iclO + T 01®y01 + - 3 - °xLl + ~T- °yll = — 


\ 


T 11 . \ X 2 . /. 

• T *yllJ + -5- v H w ULl°iOO + %00 - 


“*10 - °ylO . “rOl " V)1 'ill + a yll . T *rll\ 1 4 2 ^ i 2 ^ ° xH> 2 + °y!Q 2 ^ °*P1 2 + a jCO- Z ^ 

g + 2 S 2 I ~ SB \ 2 2 


^ * ‘Vn + U EU 

k 12 / B\ 


"^ 2 + <M> 2 1 it/- . qsio^io , ff i01 g y01 f ^cll°yll { a jtb a tf> 


^|OxOO<ryOO * 


IS 




*56.) - JL f^Xj^ |«(^0 ♦ ScIO - * ♦ W cce 2, ♦ 

, cos 2| cos 2q + a^j2£ cos { cos nj + f (^00 + ID 006 + ^01 008 + °jU. 003 ^4 co « Q + 

0^25 cos S cos n'j + fe^Tj-yu sin 2| sin 2q + 1^25 sin | sin r|) 


44 dq S& 


> - Peb 


( 21 ) 


The symbols g, tj, and £ represent the dimensionless variables in 
the x-, y-, and z-directions, and the functions $, J, and 0 are 
given as follows: 
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4 = 2cr x00 “ a yOO + ( 2a xl0 " °ylo) cos 2 £ + ( 2a x01 " a yOl) cos ^ + 
^ 2o xll “ °yll) GOS 2 | cos 2 t| + ^ 2 o xl 3 “ 0 yb^ 2 £ cos | cos t\ 

1 = 2ayoo - a x0 o + ^2a yl0 - Oxlo)c° s 2 £ + ^ 20 yOi “ °x 0 lj cos 2 r\ + 
j^20yii " °xll) cos 2 t cos 2 t) + ^2ayb “ °xb) 2 £ cos £ cos H 

0 = T-yy i ■] sin 2£ sin 2r\ + T X y^ ) 2t, sin g sin T} 


> ( 22 ) 


The conditions which the time derivatives of the coefficients must 
satisfy axe: 


an 

_ dn . 

_ an 

_ an 

an 

_ an 

an 

an 

an 

Aoo 

*°xlO 

^xOl 

S6 xll 

<*xb 

^ & y00 

^yio 

S6 yOl 

S& yil 


_ an _ 

an 

an 

_ an _ 

an _ 

an _ . 

= dn 

an 



^xyll 

^xyb 

ae 

1 

0 

£ 

au-n af 

av 0 i 

avu 


= “ = 0 (23) 

AlI 

When the operations indicated in equations (23) are performed, 
the following system of 19 ordinary differential equations is obtained: 

2 

-e + ~ v u ¥ n - I x0 o = 0 (24a) 

*“10 “ w lAi + " ^10 = 0 (24b) 
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3t ■ Jt 2 . Oxll 4Cfyll 

p U H ' 16 w n w n - 1 — + T? = 


4e k& 


ii 2 h. • 0x5 fj.0y5 

— 7 - ‘W-I -1 — + - — X-yV, = 0 

2k b 11 12E 12E xb 


f + — Vt -,V- 


- “li^u - iyoo = 0 


xt- • °ylO 

T w u"n - IS" - Vxo - 0 


ItV, 


01 


Jtf 

8 


* 11*11 + 


1*0x01 

2E 


“ %01 - 0 


it 


Oyll _ ^OjcH 


2 V H ' 16 " ~W + 1 i“ " = 0 


it 2 h . °yt) t*°xb 

2k b W H - S HT " ^ = 0 


Jt/ • , • \ 3t*^ • 1 + U • 

- f(**ii + * 11 ) + -g- W 11 - - Vn = 0 


3t 2 h ■ , 1 + Li • 

V,, + — = T 


12 b 


’ll 


6e 


T xyb + I xyb 0 


('24d) 

(2^) 

(24f) 

(2%) 

(24h) 

(sin) 

(24j) 

(2te) 

(sin) 


°xll = °yll = Txyll 

(2ifm) 

hbo xOO = 

(2ita) 

°y00 = Oxio = Oyoi = 0 

(Site) 
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q xll + 
2 


OylO 


fzii + 
2 



*11^x00 


°x01 

2 


q xii + 

4 


°ylO _ °yll + T xyll\ 
2 " 4 2 J 


+ 


1 

6 



+ a, 


yb 



= 0 


(24p) 


The integrals I in equations (24) are for various combinations of the 
stress parameters and are defined as follows (note that these eqs. do 
not contain any rates): 


1x00 = 

4J 

r 

F® dV 


3« J 

V 

^lO = 

TS J 

P F® cos 2fj dV 

TT 


3* 

V 

1x01 = 

■hj 

r F® cos 2q dV 

V 


3 * 

V 

^Kll = 

1 

2 i 

f F® cos 2g C06 2 t) dV 


3* J 

V 

•^xb ” 

-J 

F®2£ cos | cos T] dV 


3it 2 ^ 

T 

•Wll 

= 1 
0 

r 6F8 sin 2| sin 2i} dV 


3* 

0 v 




Syoo 

-ft 

II 

" Ff dV 
V 

^ylO 


n 

Ff cos 2£ dV 
V 

ZyOl 

I CVJ 

H & 

II 

r 

FJ cos 2i] dV 
V 


-ft 

11 

P 

F$ cos 2g cos 2ij dV 
V 


r 

3ji 2 A 

Ff2£ cos £ cos T] dV 
r 

I xyb 

1 w 

H | & 
11 

r 12F9t, sin | sin t) dV 
V 


(25) 


where dV = d| dTj d£, and the integration is to he performed throughout 
the volume of the plate - that is, over £ from -l/2 to l/2 and over f- 
and T) from -rt/2 to fl/2. 

The system of equations (eqs. (24)) contains all the -unknown coef- 
ficients in the stress-rate and displacement -rate distributions . These 
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equations conceivably can be solved for all these coefficients. In this 
investigation, however, equations (24b), (24f), and (24h) are not used 
since the desired information on lateral deflection, unit shortening, 
and stress distribution can be obtained without them. 

In order to start the numerical solution, the initial conditions 
or the conditions at time zero Just subsequent to loading and before 
creep has commenced are required. These conditions are discussed briefly 
in the next section. 


Initial Conditions 

Because of the nature of the load-application process assumed in 
this investigation, the initial conditions for the creep solution - that 
is, the conditions just subsequent to loading and before creep commences - 
can be obtained from a solution of the problem in which creep is neglected. 
If the applied load is below the elastic buckling load of the plate, it is 
necessary to take into account an initial imperfection of the plate in 
order to initiate lateral deflection due to creep. If the applied load 
is above the elastic buckling load, it may be possible to disregard 
initial imperfections. In this report, however, an initial imperfection 
is assumed to exist in either case. 

The Initial conditions are derived in appendix A. The results of 
the derivation are contained in equations (A7n) and (a 8) to (All) . 


Preparation for Numerical Solution of Equations 


Simplification of differential equations .- Certain of equations (24) 
are now solved to obtain w^i and the stress -rate coefficients in terms 

of the stress coefficients. Equations (24c) and (24g) yield 


*xQl = w lAl " 2I x0lj 

°yio = E (^f V11 “ ^10 j 


j 


(2 6 ) 
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From equations (24e) and (24j) 


i _ Eit 2 h „ 
°xb ~ w -> 


• 12E / T i> \ 

^11 - + Vlytj 


2(1 - n) b 




^ ■ 2 (f- „> £ i ' 11 ■ + “ Ixb ) 

From equation (24z) 


E h 


r xyb ' 1 + J 2 b + 


(27) 


( 28 ) 


The initial conditions give cr^-^ = = T xyll 1111116 zero (see 

eq. (A7n) ) , and equation (24m) states that cr x ]j_ = Oyll = Txyll* 
Consequently, 


CT xll ~ °yll “ T xyH 


(29) 


for all time. When equations ( 29 ) are considered, equations (24d) , 
(24i), and (24k) can be solved to obtain: 


^xyll “ E ( I xll + ^yll + I xyll) ^ 

Equations ( 26 ), ( 27 ), ( 28 ), and (30) are nov used in equation (24p) to 
find v-q . The result is: 


^n^oi + ^io) + 2 b 1 - + I yb) ~ 1 + (i Ix yb_ ^ 

jt 2 P , °x00 , a x01 + a ylO , a cr 
— V 11 + — + 2E + “ 


The history of lateral deflection and stress distribution can be 
obtained by solving simultaneously equations ( 3 l)> ( 26 ), ( 27 ), ( 28 ), 
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and (30) . Unit shortening may he obtained by utilizing equation (24a) . 
The initial conditions are given in appendix A. The solution of this- 
system requires considerable numerical calculation since, in general, 
the required integrals (eqs. ( 25 )) have to be evaluated numerically. 

For two cases, however, it is feasible to carry out the integrations of 
equations ( 25 ) analytically; thus the calculation is simplified. These 
two cases, represented by n = 0 and n = 1 in the creep law (see 
eq. (11)), are considered in this report. The case of n = 0 is equiv- 
alent to a material having an exponent of one in the uniaxial power creep 
law, that is a linear viscoelastic material, and is comparable to the 
solution in reference 3* The case of n = 1 is equivalent to an expo- 
nent of 3 in the uniaxial power creep law and, of course, represents a 
nonlinear material. 


Dimensionless form of equations .- It is convenient to write the 
system of differential equations and initial conditions in dimensionless 
form. For example, in equation (31) the numerator and denominator on 
the right-hand side may be divided by e cr , and both sides of the equa- 

AECT cr 2r Hct k_1 

tion may be divided by . The remaining equations of interest 

3 n 

may be nondimens ionalized, and the resulting system of differential equa- 
tions is: 


aw 

dr 


W ( I x01 + !ylo) + 2 

JL - n( Ixb + ^ 1 + n lx &_ 

3(l - ^ 2 ) T .f 
4 

- , , °x01 + *ylO 

+ °x00 + 1 + 2 


f^oi . 5(1 - „ ag . 2lx01 

df 4 dT 


da. 


r yio = 5 (i - h 2 ) w aw _ 2 - 


dT 


dT 


2I yio 


dT xyb _ 3(1 - n) dW 6 - 

dT ~ 2 dT ~ 1 + (a Txyh 


( 32 a) 

( 32 b) 

( 32 c) 

( 32 d) 


do xb 3(l + M-) dW 

dT 


12 (~ 


1 - 


^Ixb 


+ M-Iyb) 


( 32 e) 


dT 


2 
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&T 


3(1 + d) dW 12 /- 


dw 12 /- - \ 

3 ; ■ rr^v ^ + > 


dr 


'WM. = _/V 


dr 


^xll + %•!! + Ixyll 


) 


(32f) 

(32g) 


In equations (32); 



T = the dimensionless stress 

3 n 


coefficients (represented by the symbols with bars) are the stress coef- 
ficients divided by 0 cr , and the dimensionless integrals I are the 


2nv+k-l 


integrals I divided by 


e cr XEcT cr kt ' 




The initial conditions (eqs. (A7n) and (a 8) to (All)) can be written 
in dimensionless form as follows: 


T xyll - 0 

Vi - Vo - 5(1 ~ “^ (w 2 - W t 2 ) 


V - V * - 1 - !— } -(W - Ml) 

V* - - 3<1 g - ^ (W - Mi) 


W 5 + W 


8 


3(1 - n 2 ) 


(1 - P) - W ± £ 


8 


3(1 - H 2 ) 


Wi = 0 


J 


(33) 


Numerical methods are available for solving systems such as equations (32). 
Two of these methods; the modified Euler method and the Runge-Kutta method; 
are described in reference 11, and the latter was used to obtain the cal- 
culated results presented in this paper. The pertinent integrals In 
equations (25) have been evaluated for n = 0 and n = 1 and the results 
are presented in appendix B . 
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Unit Shortening and Effective Width 

One interesting feature of this analysis is the possibility of cal- 
culating the history of the unit shortening of the plate. By combining 
equations (24a) and (A7a) and writing in dimensionless form, the following 
expression is obtained for e, the ratio of the unit shortening at a given 
value of r to the critical unit shortening ' for the plate: 

e = 5 ( 1 - ^ (w 2 - W ± 2 ) + P - 1^00 dr (34) 

When the lateral deflection coefficients and stress coefficients are 
known from the solution of equations ( 32 ), the history of e can be 
calculated from equation (34). 

Closely associated with unit shortening is the concept of effective 
width. An effective width for a plate subject to creep may be defined 
analogous to the effective width of a buckled elastic or elastic-plastic 
plate. That is, the effective width at a given time and unit shortening 
is that width of plate which would support the actual load on the plate 
if the stress were uniform and equal to the stress obtained from the 
creep law at the given time and at a strain equal to the given unit 
shortening. In symbols this definition can be written as: 


(35) 

b "cl 

where b e is the effective width of the plate, and a c j Is the stress 
obtained from the creep law at the unit shortening and time associated 
with the actual plate. The proper value of a c j In equation (35) can 

be found by substituting the appropriate values of e and t from equa- 
tion (34) Into the uniaxial creep law, which, in dimensionless form, is 


e = 



(36) 


Equation ( 36 ) has only one real root for cr c j, and this root is negative 

no matter what integer value n takes. This root of equation ( 36 ), 
then, is the required value of a c j In equation (35) for a given com- 
bination of t and e. 
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RESULTS AND DISCUSSION 


Brief calculations were made of the lateral deflection history at 
the center of the plate as a function of the time parameter t for a 
plate composed of a linear viscoelastic material. For a plate with a 
cubic uniaxial creep law, more extensive calculations were made of his- 
tories of lateral deflection, unit shortening, and effective width as 
functions of r. These calculations consisted of solving numerically 
for the lateral deflection and stresses from equations (32) and initial 
conditions (eqs. (33)) • For most of the calculations involving a non- 
linear material, the lateral deflection and stresses were calculated by 
using the Runge-Kutta method and an IBM type 650 electronic data proc- 
essing machine. Then it was possible to calculate unit shortening from 
equation (3^) and effective width from equations (35) and (36 ). 

In the calculations for a nonlinear material the applied load yas 
varied from 0.4 to 1.2 times the elastic buckling load. An initial 
imperfection of 0.01 times the plate thickness was assumed for each 
value of applied load. In addition, for an applied load of 0.8 times 
the buckling load, three other initial imperf ect ions were assumed, ranging 
from 0.001 to 0.03 times the plate thickness. Some discussion of these 
results is presented in the succeeding paragraphs. 


Lateral Deflection 

The results of the lateral deflection calculations are shown in fig- 
ures 2 and 3 where lateral deflection at the center of the plate divided 
by plate thickness is plotted against the time parameter t. In figure 2 
a comparison is made between the results of the present calculations and 
the theory of reference 3 for a linear viscoelastic material (n = 0) . 
These calculations were made for an applied load equal to 0.4 times the 
elastic buckling load and an initial imperfection at the center of the 
plate of 0.01 times the plate thickness. The curve labeled "Present 
theory - small-deflection" was calculated with the large-deflection terms 
in equations (8) neglected. This result should coup are directly with 
that of* reference 3 which also is a small-deflection analysis. The dis- 
crepancy can be attributed to the fact that the uniaxial and polyaxial 
creep laws used in reference 3 are compatible only for an incompressible 
material and, therefore, the application of .the theory of reference 3 
to a compressible material is not strictly correct. The results in fig- 
ure 2. are typical . and give a qualitative indication of the comparison 
for other applied loadings and Initial imperfections. Finally, figure 2 
illustrates the divergence of the small -deflect ion theory from the large- 
deflection theory. 
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Results for a plate material having a cubic uniaxial creep law 
(n = l) are shown in figure 3 . A comparison between small- and large- 
deflection theory for the nonlinear material is also shown in figure 3 * 
Various applied loadings and, for one loading, various initial imperfec- 
tions are included in this figure. 

Results are plotted up to a lateral deflection of twice the plate 
thickness, although calculations were carried considerably beyond this 
point in some cases. Beyond a lateral deflection of twice the plate 
thickness, the accuracy of the assumed deflection shape in approximating 
the deflection at the center of the plate becomes questionable. 

In theoretical analyses of the creep behavior of columns, it has 
been the usual practice to define a collapse time as being that value 
of time at which the lateral deflection or lateral-deflection rate becomes 
infinite. A column made of a linear viscoelastic material does not exhibit 
such a collapse time but a column made of a nonlinear material does. Sim- 
ilarly, the results of the plate calculations for a linear material do not 
yield a collapse time. The lateral deflection curves shown in figure 2 
continue to increase with t but become infinite only for infinite t. 

The plate calculations for a nonlinear material do yield a collapse time 
if small-deflection theory is used. The curve for small -deflection theory 
shown in figure 3 illustrates this situation. This curve approaches 
infinity asymptotically as indicated by the vertical dash-dot line. This 
collapse time does not occur, however, until long after the small- 
deflection theory has become invalid. The method suggested in refer- 
ence 5 also yields a collapse time. The calculations from the present 
theory for a nonlinear material and based on large-deflection theory, on 
the other hand, do not yield a finite collapse time. The lateral deflec- 
tion in this case increases with time, but infinite deflection does not 
occur until infinite time. 

__ The results shown in figure 3 for various initial imperfections at 
P = 0.8 indicate that initial imperfections may have a significant 
influence on the plate creep deflections. The influence of applied load 
on creep deflections, however, is considerably more important than that 
of initial imperfections. Results similar to these have been reported 
in column creep studies (ref. 12 ). 

The curves in figure 3 for various initial imperfections at P = 0.8 
are very nearly parallel. That is, the curves for the larger values 
of can be obtained closely (not exactly) by a simple translation to 

the left. At a given value of W the slopes of the curves for W 
plotted against r for a given value of P appear to be practically 
independent of Wq. In the lower part of figure 3 the t scale is 
magnified to illustrate this situation. 
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Unit Shortening and Effective Width 

The calculated unit shortening of the plate as a function of t 
is presented in figure 4. Note that the t scale is magnified for 
clarity in the lower part of this figure. The dashed lines were computed 
from the uniaxial creep law which can he written in dimensionless form 
as: 


6 = "( a + I ~°J (37) 

If the initial imperfection is sufficiently small and the applied load 
is less them the elastic buckling load, the unit shortening curves start 
out very nearly identical to the creep law and diverge smoothly from the 
creep law. The unit shortening curves become identical to the creep 
curves if the initial imperfection. goes to zero. Thus when t is small, 
the creep law can provide a first approximation to the unit shortening 
as might be expected. For applied loads higher than the elastic buckling 
load, such an approximation may not be valid as the unit shortening curves 
start out above the creep curves . 

A plot of effective width divided by actual width of plate against r 
is presented in figure 5* Again, in the lower portion of the figure the t 
scale is magnified for clarity. It is seen that creep may cause a signif- 
icant decrease in the effective width of a plate. 


Stress Distribution 

In figure 6 are presented the middle-surface stresses along the 
edges of a plate with P = 0.8 and Wi = 0.01 at two values of t. 

As a result of creep there can be considerable redistribution of middle 
surface stress as time increases. Along the edges of the plate which 
are subjected to the applied load there is a tendency for the stresses 
at the edges to increase and the stress in the center to decrease. Fur- 
thermore, as creep progresses, significant stresses can arise along the 
so-called unloaded edges of the plate provided these edges are constrained 
to remain straight. The middle -surface stress distribution in a square 
plate undergoing creep, then, resembles that in a buckled elastic plate. 


Collapse Time 

As previously mentioned, the large -deflect ion analysis in this report 
does not yield a finite collapse time for the plate of a nonlinear material. 
It is known from experiment (ref. l) that plates supported in V-groove fix- 
tures do. collapse in finite time when subject to creep. The V-groove 
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fixtures approximate the simply supported edge condition but do not pro- 
vide the straight edge constraint assumed for this analysis. This assump- 
tion of straight edges seems to be more realistic than the V-groove sup- 
port in approximating a portion of sheet lying between stiffeners in 
certain multiple bay stiffened panels, for example. Although a plate 
which satisfies the edge conditions assumed in this analysis may last 
longer than a V-groove supported plate, it does not seem reasonable that 
such a plate would never collapse. These matters probably cannot be 
clarified without further experimental work. 

An actual plate or col umn , of course, always has finite deflections 
and strains even after collapse. The usual definition of collapse time 
appears reasonable for a column because calculations show that infinite 
deflection rates or strain rates are approached while the theory is still 
supposed to be a good approximation of reality. The results of the pres- 
ent investigation indicate that this situation is not valid for plates. 
Perhaps same additional considerations are required before a satisfactory 
theoretical description of plate creep collapse can be given. 

There is a similarity between the results of the present analysis 
and the results of reference 10. In that investigation an analysis was 
made of a plate which was buckled in the elastic range and then compres- 
sed into the plastic range. It was found that no maximum load resulted 
from the analysis. In reference 10 possible changes in buckle pattern 
which occur in experiment were not taken into account. Such changes 
have been ruled out of the present analysis, and this assumption may be 
satisfactory for a square plate. Some consideration probably should be 
given to this phenomenon, however, in the study of the creep of rectangular 
plates with aspect ratios other than unity. Another possible refinement 
in the present analysis is the use of an improved stress distribution 
through the thickness of the plate. In addition, modification of the 
creep law or the values of the constants in the creep law may be required 
at large creep strains. 

Even though a theory does not yield quantitatively useful results, 
it may aid in finding parameters which correlate test data. The present 
theory suggests that a useful correlation might result from a plot of P 
against' the time parameter t with the experimental collapse time sub- 
stituted in place of t. In such a procedure perhaps It could be assumed 
that Initial imperfections would not be widely different for flat plates 
manufactured by similar methods. It Is interesting to note that r Is 
a simple parameter which includes material properties and plate geometry 
and can be obtained directly from the creep law without recourse to an 
elaborate theory. 
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CONCLUDING REMARKS 


An analysis has been made of a square plate composed of a material 
subject to creep. The plate is assumed to have small initial imperfec- 
tions. The material is assumed to follow a generalized power creep law 
which accounts for the biaxial state of stress. Calculations have been 
made for a value of the exponent in the generalized creep law which cor- 
responds to an exponent of 3 in the uniaxial power creep law and for a 
linear viscoelastic plate. 

t 

Calculations based on small-deflection theory and a nonlinear material 
yield a collapse time - that is, a finite time at which the lateral deflec- 
tion becomes infinite. The more refined calculations based on large- 
deflection plate theory, however, do not yield a finite collapse time. 

Test results show that plates supported in V-groove fixtures do collapse 
in finite time. The assumptions underlying the present theory, however, 
seem to be more realistic than V-groove supports for certain types of 
plates which might be encountered in practice. Thus, the theory indicates 
the possibility that certain practical plates may carry a given load sig- 
nificantly longer than geometrically similar V-groove supported plates. 

It is found that although initial imperfections have a significant 
influence on creep deflections, applied load, has a much stronger influ- 
ence. Creep in a plate can cause considerable redistribution of the 
middle-surface stresses. Along the loaded edges of the plate the stresses 
tend to increase at the ends and decrease in the center. In addition, 
significant stresses can grow along the unloaded edges provided these 
edges remain straight. Finally, as creep progresses, the effective width 
of the plate Is reduced. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 23, 1958. 
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APPENDIX A 

DERIVATION OF INITIAL CONDITIONS 


In this appendix the initial conditions for the creep solution are 
derived. These conditions define the state of stress and strain in the 
plate just subsequent to the load application but prior to the beg inning 
of creep. Thus, it is an elastic solution with no creep which is sought 
here. 


Reissner's Principle 


An approximate elastic solution can be obtained by the use of 
Reissner's variational theorem. (See ref. J.) The form of the theorem 
appropriate for the present problem can be stated as follows: The state 

of stress and displacement which exists in the plate is determined by 
the variational equation SUg = 0 where 


% 


ph/2 

ob/2 

ob/2 

> -h/2 ^ 

-b/2^ 

-b/2 







+ 


2 / \ ‘ 
— + | l( ^i \ 

dy 2 \byJ 2\dy J 



*2 
o w-l 


dy 2 J 


°y + 



Sv + 5w dw 
chc chc Sy 



In equation (Al), 


w^ represents the initial imperfection in the plate. 


A Rayleigh-Ritz procedure is used to obtain a solution. Spatial 
distributions^ which contain arbitrary coefficients, are assumed for the 
displacements and stresses in the plate. These expressions are substituted 
into the equation for Hr, and the variation of IIr is equated to zero. 
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This procedure is the sane as equating to zero the partial derivatives 
of % -with respect to each unknown coefficient. There results a system 

of s imu ltaneous linear algebraic equations which determine the u n k n owns. 


Ass ume d Displacements and Stresses 
The displacement and stress distributions are taken as follows: 


, l 2ity| . 2rtx 

u = -ex + bluj_o + u n cos -raisin 


v = fy + b(v Q1 + v xl cos ^*)sin ^ 


. jtx ity 

w = bwTi cos — cos 

J - L b b 


i> 


(A2) 


°x “ °x00 + °xlO 008 + °x01 C08 + °y 1 1 cos cos ^ + Ojcb ^ cos ^ cos 22 

a y = Oy00 + °ylO = 08 + CyOl cob ^ + Oy U cos ^ cos ^2 + ayb cos ^ cos S l (A3) 


T — = T xyli 


xy 




These distributions are the same as equations (19) and (20) except the 
dots are omitted. Here, of course, the assumption of a linear variation 
■ of stresses in the z-direction is consistent with elastic material 
behayior. The initial imperfection is assumed to be in the same shape 
as 'the lateral deflection in equations (A2), that is: 


, tcx icy 

wj_ = b v-n-t cos — cos — 
b b 


(A k) 


Algebraic Equations for Coefficients 

When equations (A2), (A3), and (A4) are substituted into equation (Al) 
for II R and the indicated integrations performed, there is obtained: 



NACA TN ^398 


29 


Hr = lab 4 


-°xOO e + °yCO f + y( V ll 2 * v lli 2 )^< 


„ , °xlO " °ylO 

a xOO + CT yOO p + 


°x01 ~ °y01 0 xll + °yll T xyll^ 


■ 2 


+ * °xl0 u 10 + °yOl v 01 + 


g xii u n + p yii v n u n + v n 


°yb - 2r xyb 


r xyllj + ^ ^( v ll “ w lli) ( p xb + 

) - i(«xoo 2 + °y 00 2 ♦ S3 S- l ^ * 5s! ; VI s 


0 xl 1 2 + OF-. ~ ^ <t.2-i_„_2 


F yll , °xb + °yb 1 , \i(_ _ , 0 xlO°ylO + p x01°y01 . 

t + £ J + £ IWyOO + 5 + 


°xll q 'yll + p xb °yb\ _ 1 + n 

)l to/ t ottt 


12 


12E r T xyll + T xyb £ 


- Peb 


(A5) 


The conditions which the unknown coefficients must satisfy are: 


Se 


an E 

ai E 

an E 

3% 

3% 

3n E 

^°x00 

^ a xl0 

Sa x01 

Sa xll 

So xb 

^°yOO 


... ^R 

8n R 

_ Sn R 


8n R 


aa y01 

^°yll 

3lyU 

^xyll 

^xyb 


Sn R 8n R 8 it r 3n R Sn R Sii] 


L R 


Ho ^11 Sf ^01 ^11 


= 0 


(AS) 


After the operations indicated in equations (A6) are performed, some 
algebraic manipulation leads to the following 19 equations: 


+ ir( v n 2 ” v lli 2 ) ■ 


p x00 

E 


= 0 


(A7a) 
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50110 “ 16 ^ ll2 " Vlli2 ) + Ie - ~ 0 

(A?b) 


jt 2 / p 2 \ ff x oi n 

Jg("rL 2 * w ni 2 ) 2E ---° 

(A7c) 


| U 11 - ^(“ 11 2 - w iii 2 ) - 0 

(A7d) 


it 2 h/ . \ °xb " ^ a yb _ A 

2k b\ 11 111 ) 12E 

(A7e) 


f + ^ . W ) + ^ . 0 

(A7f) 


s(^ 2 ' ^ " 0 

(A7g) 


™01 - fg^ll 2 - w lli 2 ) + 2 E - 0 

(A7b) 


| V 11 " ^2 ( w ll 2 ” w lli 2 ) = 0 

(A7i) 


it 2 h/ \ a yb “ ^ _ n 

i4 b V 11 “ I 2 E 

(A7J) 


£ ftii - 'in) + « “ ^ ■ ° 

(A7k) 

^xOO 

+ .rf 1 + ^f 2 ) + 1 f (°*b + «*> - 2 T *Jf) = 0 

(A7D 


^xOO = " P 

(A7m) 

0 yoo 

= °xlO = 0 yOl = a xll = °yll = T xyll = 0 

(A7n) 



NACA IN 4398 


31 


Solution of Algebraic Equations 

Not all of equations (A 7 ) are required to obtain the information 
needed in this investigation. The procedure is to solve for the stress 
coefficients in terms of w-^ and w-m and then to substitute these 

expressions into equation (AJl ) . Thus, there is obtained an equation 
from -which Vjjl at time zero can be calculated. 

From equations (A7c) and (A7g): 

°xoi = °yio = ^( v n 2 " w ni 2 ) ^ A8) 

The s imult aneous solution of equations (A7e) and (A7j ) yields: 

“*> ‘ ‘ 2(1”! n) I ( Wl1 " Wlli ) (A9) 


From equation (A7k) the following is obtained: 


T xyb = 




(A10) 


When equations (A 8 ), (A 9 ), and (A10) are substituted into equation (A7Z) 

p 

and use is made of the fact that o^oo = (eq. (A7m)), the result 

* bh 

can be written as follows: 


"n 3 + "‘111 



8 


3(1 - 



ill 


= 0 (All) 


where P is the ratio of the applied load to the elastic buckling load. 
The positive real root of equation (All), which can be obtained by trial 
and error, gives the initial condition for - that is, the value 

of w-j^ immediately subsequent to the load application. The initial 

conditions for the stress coefficients can then be calculated from equa- 
tions (A 8 ), (A9), and (A10). 
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APPENDIX B 

RESULTS OF EVALUATION OF INTEGRALS 
REQUIRED FOR CALCULATIONS 


In order to perform the calculations for this report the integrals I 
in equation (32) and equation (3^0 must he evaluated for the proper values 
of n, the exponent in equation (ll) for the function F. The results 
are as follows: 

For n = 0 (linear viscoelastic material) 
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For n = 1 (cubic uniaxial creep law) 
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Figure !•- Plate geometry and coordinate system considered in analysis. 
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Figure 2.- Comparison ' of histories of lateral deflection at center of 
plate calculated by three theories for a linear viscoelastic mate- 
rial. P = O.if-j = 0.01; n = 0; n = 0^3- 
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Figure 3>- History of lateral deflection at the center of a square 
plate made of a nonlinear material, n = 1; n = 0.3* 
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Figure 5- - History of effective width of square plate made of a non- 
linear material, n = 1; |i = 0-3- 





